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We consider a scalar charge travelling in a curved background spacetime. We calculate the quasi- 
local contribution to the scalar self-force experienced by such a particle following a geodesic in a 
general spacetime. We also show that if we assume a massless field and a vacuum background 
spacetime, the expression for the self-force simplifies significantly. We consider some specific cases 
whose gravitational analog are of immediate physical interest for the calculation of radiation reaction 
corrected orbits of binary black hole systems. These systems are expected to be detectable by the 
LISA space based gravitational wave observatory. We also investigate how alternate techniques may 
be employed in some specific cases and use these as a check on our own results. 



I. INTRODUCTION 



There has been considerable recent interest in the modelling of astronomical and cosmological gravitational wave 
sources. This interest has been spurred on by the recent and upcoming construction of both ground and space based 
l/^ ' gravitational wave observatories. Ground based detectors such as LIGO VIRGO and GEO600 H have already 
reached the data collection phase. The primary space based gravitational wave observatory, LISA ^ - the Laser 
Interferometer Space Antenna - is still in the design and planning phase. A joint NASA/ESA mission, it is expected 
^h' to launch within the next few years. It will detect gravitational waves in the frequency band 0.1-O.OOOlHz, a band in 
^ . which a considerable number of interesting gravitational wave sources are expected to be found. 

In order to extract the maximum amount of information from these gravitational wave detectors, data analysis 
techniques such as matched filtering may be employed. To improve the effectiveness of this data analysis, it is 
necessary to have accurate predictions of the waveforms of the gravitational radiation emitted by candidate sources. 
One of the primary candidate sources of such a gravitational wave signal falls under the category of large mass binary 
' systems. In particular, extreme mass ratio systems such as that of a compact solar mass object (mass m) inspiralling 
\^ \ into a black hole of lO'^ to 10* solar masses (mass M) are expected to be detectable by LISA out to Gpc distances. 

■ As a prerequisite to the calculation of the waveform templates for such sources, it is necessary to model the orbital 
O^l ' evolution of the system. In the case of extreme mass ratio systems, we can do so by making the approximation that 
the smaller mass is travelling in the background spacetime of the larger mass. However, the compact nature of the 
smaller mass means that it itself will distort the curvature of the spacetime in which it is moving to a non-negligible 
, extent. As a result, it does not follow a geodesic of the background spacetime of the larger mass, but rather, it follows 
' a geodesic of the total spacetime of both masses. In this sense the motion of the smaller mass is affected by the 
presence of its own mass, i.e. the smaller mass is seen to exert a force on itself. This is called the self-force. Clearly, 
the computation of this self-force is of fundamental importance to the accurate calculation of the orbital evolution of 
such binary systems and hence to the prediction of the gravitational radiation waveform. 

Because of this non-linear behavior, it is extremely difficult, if not impossible, to exactly solve the equations of 
motion of such a system. However, due to the extreme mass ratio involved, the deviation of the object's trajectory 
from the background geodesic is small over sufficiently small time-scales (less than the natural length-scale of the 
background). In this case, to linear order in m, we can model the deviation from the background as a spin-2 field 
generated by the particle and living on the background spacetime. This field is then seen to couple to the smaller 
object to cause the self-force. 

In the present work, we will focus not on the gravitational sclf-forcc described thus far, but rather on the analogous 
but simpler scalar self-force. Instead of considering the smaller mass, m, to be generating a gravitational field, we 
take it to have a scalar charge, q. This charge couples to a massless scalar field which is then the cause of the self- 
force. Otherwise, we leave the problem exactly as posed above. In making this change, the resulting analysis remains 
extremely similar to the gravitational case, but the algebraic details of the calculations become simpler. This will 
allow us to explore the problem and develop techniques which we will at a later date apply to the physically more 
interesting gravitational case. 



* Electronic address: 'adrian.ottewiU@ucd.ie] 
t Electronic address: barry.wardellOucd.ie, 



2 



An expression for the electromagnetic self-force was originally obtained by DeWitt and Brehme 0] (later corrected 
by Hobbs Q) and more recently recovered by Quinn and Wald Q- In addition, Ref. Q derives an expression for the 
gravitational self-force, as was done previously by Mino, Sasaki and Tanaka Using a generalization of his axiomatic 
approach, Quinn also obtained an expression for the scalar self-force 0]. These results show how the self- force may be 
expressed in terms of an integral of the retarded Green's function over the entire past world-line of the particle. This 
retarded Green's function possesses two terms, one representing geodesic propagation, the other a tail-term arising 
from back-scattering by the spacetimc geometry. (Solutions to wave equations on most curved manifolds depend not 
only on Cauchy data directly intersecting the past light cone, but also on Cauchy data interior to that intersection. 
The portion of the field that propagates in the null directions is called the direct part, while the portion of the field 
that propagates in the timelike directions inside the light cone is called the tail.) Clearly, the tail part of a particle's 
field can interact with the particle, leading to a contribution to the self-force. 

Although in principle this formal solution gives the desired result, in practice the calculation of the Green's function 
poses a formidable challenge. There has been a significant amount of work in this area in recent years. Various 
techniques have been emp loyed. There has been considerable recent success using mode-sum approaches, amongst 
others [13, El, [13, [3, [3- Other work has focused on a matched expansion approach and also on obtaining exact 
solutions under specific circumstances. For the Schwarzschild background, Wisem an jl5| found an exact solution for 
the scalar self-force on a particle held at rest. Anderson, Hu and Eftekharzadeh [l^ built on this work along with 
other results to derive an expression for the self-force on a particle held at rest, then released and allowed to fall 
radially inward. Anderson and Wiseman [l7| also examined the matched expansion approach under the two cases 
of static particles and circular geodesic paths in a Schwarzschild background. For the case of a gravitational field, 
Anderson, Flanagan and Ottewill [l^ have evaluated the quasi-local contribution to the self-force up to O (Ar^^). 
Comprehensive reviews of the radiation-reaction and self-force problem are given by Poisson [Toj and by Detweiler 

In this paper, we use the results of Dccanini and Folacci [2^ as a basis for our own work. In section |TT] we first 
recap the work of Quinn [^ to define what we mean by the quasi-local part of the self-force and then we derive an 
expression for it and for the equations of motion in terms of the Green's function. Next, in section [TTIl we use this 
result to calculate an expression for the quasi-local portion of the scalar self-force in terms of the expansion coefficients 
given by Ref. ^2], and in terms of the proper time to the matching point At and the 4- velocity of the particle, u". 

By making the assumption of motion in a massless field on a vacuum background spacetimc, we show in section HVl 
that the resulting self-force expression simplifies considerably and in fact vanishes up to and including third order in 
At. 

In section|Vl we examine some particularly interesting and simple cases in both Schwarzschild and Kerr spacetimes. 
We then show that using a corrected form of the results of Ref. [l3|, we obtain identical results, order by order for 
these special cases. This confirms the consistency of our results and the validity of our approach. We also present 
corrected results for the higher order terms in Ref. [l7| and give updated truncation error graphs. 

Throughout this paper, we use units in which G = c = 1 and adopt the sign conventions of [1^. We denote 
symmetrization of indices using brackets (e.g. (a/3)) and exclude indices from symmctrization by surrounding them 
by vertical bars (e.g. (a|/?|7)). Roman letters are used for free indices and Greek letters for indices summed over all 
spacetimc dimensions. 



II. THE QUASI-LOCAL SCALAR SELF-FORCE AND EQUATIONS OF MOTION 

Consider a particle of scalar charge q and mass m travelling in a curved background spacetimc. We treat this as 
a point particle and will refer to it as a particle throughout this paper. However, the arguments presented here are 
equally valid for non-point objects provided the object's extended structure does not affect its center-of-mass motion. 

The scalar field generated by the particle propagates along null geodesies, while the particle itself will (approxi- 
mately) follow a time-like geodesic of the background. The field then interacts with the particle, effectively causing a 
force, called the scalar self-force. This self-force may be expressed as a sum of local and non-local parts [0| 

r = Q - a'"") + I {R''^up + Rpyu-'u'') + QmLd " ^ (1 " ^0 +1™ f ' VG.et (x, x') dr' 

(2.1) 

where u" is the 4- velocity of the particle, a is the 4-accelcration, a = I7 is the derivative of the 4-acceleration with 
respect to proper time and Gret is the retarded scalar Green's function. Here, we have trivially extended the expression 
given in Ref. (9| and Ref. [l9| to include terms corresponding to the field mass mfidd and coupling to the background 
scalar curvature, ^. 
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In this form, the local terms are all given explicitly, so their calculation can be carried out immediately without 
posing any real difficulty. Our task is therefore to elucidate the non-local integral term. To this end, we will work 
with an expression for the self-force which just contains the non-local integral term: 

/r — e 
V^G.et (a;, x') dr' (2.2) 
-oo 

We do so with the understanding that the local terms can easily be added back in later if necessary. 

This leaves us with an expression for the scalar self-force in terms of an integral of the gradient of the retarded 
Green's function over the entire past world line of the particle. The limiting feature at the upper limit of integration 
is necessary because of the singular nature of the Green's fimction at the point x ~ x' . 

This integral would normally be evaluated using a mode-sum approach. One would calculate a sufficient number of 
the modes and sum them up to obtain the Green's function solution and hence the self-force. However, the limiting 
feature in the upper limit of integration causes great difficulties for this method. As e decreases, it is necessary to 
take an increasing number of modes to obtain an accurate solution. When e tends towards 0, the number of modes 
required becomes infinite and it is no longer possible to accurately compute the integral. However, we can avoid this 
problem if we split the integral in two parts 

/"==limgM \7''Gret{x,x')dT' +q'^ V''Gret{x,x)dT' (2.3) 

Jt-At J-oc 

where r — At is a matching point. It is chosen so that both integrals may be evaluated separately using different 
means and then combined to obtain an expression for the total self-force. The second integral in this expression 
may be evaluated, for example, by using the methods of quantum field theory in curved spacetime [l^l to perform a 
mode-sum calculation of the Green's function away from the singularity x = x' . We call the first term in Eq. (j2.3p 
the quasi-local part of the self-force. It is the computation of this integral that is the main focus of the remainder of 
this paper. 

The greatest obstacle to obtaining a solution for this quasi-local integral lies with the difficulty in the computation 
of the retarded Green's function, Gret- It is a solution of the wave equation on a curved background spacetime: 

(□, - - CR) Gret {x, x') - -4^5 {x, x') (2.4) 

In general, it proves extremely difficult, if not impossible, to find an exact solution of Eq. (|2.4p . However, provided 
X and x' are sufficiently close, we can use the Hadamard form for the retarded Green's function solution [Qj, 

Gret (x, x') ^ [x, x') {U (x, x') 5 (cr (x, x')) ~ V (x, x') 9 (-cr (x, x'))} (2.5) 

where 6- {x,x') is analogous to the Heaviside step-function (i.e. 1 when x' is in the causal past of x, otherwise), 
5 {x,x') is the standard Dirac delta function, U {x,x') and Vix,x') are symmetric bi-scalars having the benefit that 
they are regular for x' — > x and a (x,x') is the Synge 0, [111, [23 world function. For the time-like geodesies that are 
of interest to this paper, a is equal to minus one half of the squared geodesic distance between x and x': 

a(x,x') = -^s2 (2.6) 

This form of the retarded Green's function is valid provided x' lies within a normal neighborhood of x, i.e. provided 
there is a unique geodesic connecting x and .x'. Clearly, by simply choosing an appropriate value for Ar, we can 
satisfy this normal neighborhood condition and use the Hadamard form for the retarded Green's function in the ffi'st 
integral in Eq. ([Q]) .^ 

The first term (involving U{x,x')) in Eq. (|2.5[) is called the direct part of the Green's fimction and does not 
contribute to the self-force. This can be understood by the fact that the presence of S {a (x, x')) in the direct part of 
Gret means that it is only non-zero on the past light-cone of the particle. However, since the quasi-local integral in 
Eq. (|2.3p is totally internal to the past light-cone, this term does not contribute to the integral. 

The second term in Eq. (|2.5p is known as the tail part of the Green's function. It is this term that is responsible 
for the self-force and its calculation is the primary focus of this paper. 



^ Anderson and Wiseman [TtI have calculated some explicit values for the normal neighborhood boundary for a particle following a circular 
geodesic in Schwarzschild spacetime. These values place an upper bound estimation on the geodesic distance in the past at which we 
can place the matching point Ar. 
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We can now substitute Eq. (|2.5p into the quasi-local part of Eq. (|2.3p to obtain the final form of the quasi-local 
self-force: 



Jql — 1 



f V''V{x,x')dT' (2.7) 



Note that it is no longer necessary to take the limit since V {x, x') is regular everywhere. 

Unlike the cases of electromagnetic and gravitational fields, where the 4-acceleration is exactly equal to the self-force 
divided by the mass, there is one further step required in order to obtain the equations of motion in the presence of a 
scalar field. In this case, the self-force contains two components, one related to the 4-acceleration, the other related 
to the change in mass: 



r^i.^^u'^ + ma'^ (2.8) 



The 4-acceleration is always orthogonal to the 4-velocity, so in order to obtain it from the self-force, we take the 
projection of /" orthogonal to u°: 

ma" = P^f (2.9) 

where 

Pl^5i + uau' (2.10) 

is the projection orthogonal to it". 

The remaining component of the self-force (which is tangent to the 4-velocity) is simply the change in mass: 

^ = -r«a (2.11) 

III. COVARIANT CALCULATION 

Within the Hadamard approach, the symmetric bi-scalar V {x, x') is expressed in terms of an expansion in increasing 
powers of a [2^ : 

oo 

V{x,x') = J2'^n{x,x')a''{x,x') (3.1) 

n=0 

Substituting the expansion (|3.ip into Eq. (|2.4p . using the identity (J-aO'''°' — 2a and explicitly setting the coefficient of 
each manifest power of cr" equal to zero, we find (in 4-dimensional space-time) that the coefficients Vn {x,x') satisfy 
the recursion relations 

{n + 1) (2n -1- 4) Vn+i + 2 (n + 1) K-Hi;pcr^^ - 2 (n + 1) l/„+i A'^/^^^^^^a'^ + (D^ - - ^R) V„ = for n e N 

(3.2) 

along with the initial conditions 

2Vo + 2Vo.f,a'^' ~2VoA~'^/'^A^^^,^a''' + {D^^m^ -^R)Uo = (3.3a) 

Uo = (3.3b) 

The quantity A {x,x') is the Van Vleck-Morctte determinant defined as [l^l 

A (x, x') = - [-g (x)]-^/' det {-a,^,. {x, x')) [-g {x'T'^' (3.4) 

It is important to note that the Hadamard expansion is not a Taylor series and is not unique. Indeed, in Appendix [B] 
we show how, under the assumption of a massless field on a vacuum background spacetime, the terms from Vb and Vi 
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conspire to cancel to third order when combined to form V. Such cancellation is also seen, for example, in deSitter 
space-time, where for a conformally invariant theory all the Vn^s are non-zero while = 0. 

In order to solve the recursion relations (|3.2p and (|3.3p for the coefhcients ¥„ {x,x'), it is convenient to first write 
Vn {x,x') in terms of its covariant Taylor series expansion about x = x': 



Vn (X, x') = ^ —'"n{p) {X, x') 



p=0 



where the Vn{p) are bi-scalars of the form 



Vn{p) {X,x') = Vn ai...ap{x)a'°'^ {x,x') . ..a'"" {x,x') 



(3.5) 



(3.6) 



Dccanini and Folacci [g^l calculated expressions for these u„(p) up to O (c^), or equivalently O (Ar"'), where At is 
the magnitude of the proper time separation between x and x' . In Appendix El we use the symmetry of the Green's 
function in x and x' to extend their results by one order in the proper time separation, Ar, or equivalently one half 
order in a {x,x'). 

We now turn to the evaluation of the self-force. We could calculate the self-force given this form for Vn {x,x'), 
however, it proves easier to work with V {x, x') expressed in the form of a covariant Taylor expansion, i.e. an expansion 
in increasing powers of cr'°: 



V [X, X') = V [xy^ (X, X')... a^"- {X, X') 



(3.7) 



Explicit expressions for the 'yai...ap (2;), up to O (cr^^^), can then be easily obtained order by order by combining the 
relevant ai...ap as described in Appendix [Bl 

Substituting Eq. (|3.7p into Eq. (j2.7p . taking the covariant derivative and noting that the i'ai...a are symmetric 
under exchange of all indices, ai, . . . , Op, we obtain an expression for the self-force in terms of Wai...ap5 and their 
derivatives: 



fo 



QL 



-At 



Now, by definition, for a time-like geodesic. 



dr' 



a'" = -(t-t')u" 



(3.8) 



(3.9a) 
(3.9b) 



where w° is the 4- velocity of the particle. 

Also, the quantity ct'"'' can be written in terms of a covariant Taylor scries expansion (sl. [l9l. [2^ : 



•.ah ah r>a h -.a -[3 , ^ r^a h -.a -.3 



— R'' 
60 



a /3;7(5 



p b 



45 



^;a^;/3^;7^;5^0|^(^;a)5^ (3.10) 



Substituting Eqs. p.9ap . (j3.9bp and (|3.10p into Eq. (|3.8|) . the integrand becomes an expansion in powers of the 
proper time, r — r': 



/ql 



-At 



^"(r - r')° + B"(r - r')^ + C"(r - r'f 



+ D''{t - T'f + E''{t - r'f + O ((r - r')') 



dr' 



(3.11) 
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where 





= V'"- 




:a 

= Wfl ' 
P 




1 

= 


jja 


1 

= 12^' 







(3.12a) 

- W/37/"m^ (3.12b) 

+ ^I'^V^M^w'' - i^/s^i/'^w-'u^ (3.12c) 



(3.12d) 



+ \vp^sR^;^v?u'u'u'^ + Ly^^^^^-ul'u^u'u' - ^vp^secg^Vu'u'u^ (3.12e) 

The coefRcients of the powers of r — r' are all local quantities at x, but the integral is in terms of the primed 
coordinates, x' , so the coefficients may be treated as constants while performing the integration. The result is that 
we only have to perform a trivial integral of powers of t — t'. Performing the integration gives us our final result, the 
quasi-local part of the self force in terms of an expansion in powers of the matching point, At: 

/^L = -Q^ ( A"At1 + -B'^At^ + -CAt^ + -D^At^ + -E^At^ + O (At^) ) (3.13) 
y 2 3 4 5 / 

It is possible to simplify Eqs. p.l2ap - (|3.12ep further by using the results of Appendix \X\ to rewrite all the odd 
order Wai...ap terms in terms of the lower order even terms. It is then straightforward to substitute in for the Vai...ap 
to get an expression for the quasi-local self-force in terms of only the scalar charge, field mass, coupling to the scalar 
background, 4-velocity, and products of Ricmann tensor components of the background spacetime. However, we 
choose not to do so here as that would result in excessively lengthy expressions. Instead, we direct the reader to Ref. 
[2^ for a general expression for the Wnai...ap for a scalar field and to Appendix [Bl where we give the expressions for 
the Vnai...ap for massless scalar fields in vacuum spacetimes up to the orders required. 

Once an expression has been obtained for the self-force using this method, it is straightforward to use Eqs. 
p. lip to obtain the equations of motion. 



IV. SIMPLIFICATION OF THE QUASI-LOCAL SELF-FORCE FOR MASSLESS FIELDS IN VACUUM 

SPACE-TIMES 

Eq. (|3.13p is valid for any geodesic motion in any scalar field (massless or massive) in any background spacetime. 
At first it may appear that this result is quite difhcult to work with due to the length of the expressions given in Ref. 
(2^ for the Vnai...ap- However, if we make the two assumptions that: 

1. The scalar field is massless (mficid = 0) 

2. The background spacetime is a vacuum spacetime {Raf} ~ 0) 

the majority of these terms vanish and we are left with the much more manageable expressions for the Vnai...ap given 
in Appendix [Bl 

In Appendix iBl we give the expansion coefhcients w„ai...ap up to O (c^^^) for a massless field in a vacuum space- 
time. We can see immediately from these expressions that all terms involving Vafi and Vap-y identically vanish. 
That leaves expressions involving only v^p-ys and Vap-ySt- However, Eq. (jAl[) gives us a relation which allows us to 
express any odd order term (in this case, Vap-yso) in terms of all the lower order terms. This, along with the fact that 
Vai...a is totally symmetric, leads to a vastly simplified expression for the quasi-local self-force on a scalar particle 
following a geodesic in a vacuum background spacetime: 



/ql = 9' [j^vp-,s.9^'^u'<u'u^At' - - (^-«75.C;/3 - 2«;3^^,;cj 9^'Vu'u^u^At'' + O (Ar^) j (4.1) 

We see that a large number of terms in p.l3p either identically vanish or cancel each other out. In fact, at the first 
three orders in At, we get no contribution to the quasi-local self-force whatsoever. Additionally, our fourth and fifth 
order coefficients take on a much simpler form than the general case given in Eqs. (|3.12dp and ((XT^ . 
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V. SOME SPECIAL CASES 



Expression (|4.1[) is a very general result, giving the quasi- local part of the self- force on a scalar charge following 
an arbitrary time-like geodesic in any 4-dimensional vacuum spacetime. While such a general result is of tremendous 
benefit in terms of flexibility, it is also interesting to examine some specific cases whose gravitational analog are of 
immediate physical interest. We will consider two well-known black hole spacetimes, Schwarzschild and Kerr, and 
some physically probable geodesies paths. 

Although these results apply only to the non-physical scalar charge, we envisage that the techniques developed here 
could be applied to the more complicated gravitational case. This would allow us to work out the quasi-local part of 
the gravitational self-force. In conjunction with a mode-sum approximation of the far-field part of the integral, we 
would then be able to calculate the total self-force. Hence it would be possible to calculate the waveform corresponding 
to the gravitational radiation emitted from such a system. These waveforms are expected to be of great importance 
in the detection of gravitational wave sources by LISA 13, the space based gravitational wave observatory . 



A. Schwarzschild Space-time 



We begin with Schwarzschild spacetime as that is the simpler of the two cases. The line-element of this spacetime 
is given by 



ds' 



dt^ 



(5.1) 



From this, we calculate all necessary components of the Riemann tensor and its covariant derivatives and hence values 
for all the Vai...a ■ While the calculations involved pose no conceptual difficulty, there is a great deal of scope for 
numerical slips. Fortunately, they are well suited to being done using a computer algebra system such as the GRTensor 
m package for the Maple CAS 

We consider a particle with a general geodesic 4- velocity travelling in this background spacetime. At first glance, 
it may appear that there are four independent components of this 4-velocity. However, because of the spherical 
symmetry of the problem, we may arbitrarily choose the orientation of the equatorial axis, 9. By choosing this to 
line up with the 4-velocity (i.e. 9 = we can see immediately that = 0. Furthermore, the particle is following 
a time-like geodesic, so the normalization condition on the 4-velocity u^u" = — 1 holds. This allows us to eliminate 
one of the three remaining components. In our case, we choose to eliminate , which has the equation: 



r-2M 



r - 2M 



(5.2) 



We are now left with just two independent components of the 4-velocity, u'^ and u^. Substituting in to Eq. (|4.ip 
the values for the known 4-velocity components and u* along with the Riemann tensor components, we get a 
general expression for the scalar self-force in terms of the scalar charge, q, the matching point, Ar, the mass of the 
Schwarzschild black hole, M, the radial distance of the scalar charge from the black hole, r and the two independent 
components of the 4-velocity of the scalar charge, and u'^. We can then use Eqs. (|2.9p - (|2.1ip to obtain the 
quasi-local contribution to the equations of motion: 



ma, 



QL 



11200ri 



ma^L 



lOr^u'' 



r - 2M 



r - 2M 



1 + 11 {ru'f'f + 13 {ru'l')^^ + {u'-f (l + 8 {ru't'f^ 



At" 



(20r - 49M) + 8 (llr - 28M) {ru'''f + 2 (34r - 89Af ) (rw'^)'*) 
+ {u^' f ((4r - 17A/) - 8 (31r - 59M) (ru"^)^ - 30 (14r - 29M) (rw"^)^) 
-16r {u"^ f (l + 15 {ru'f'Y^ 



Ar^+OfAr^) 



O (At^) 



(5.3a) 
(5.3b) 
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3g^M^ 
11200rio 



ma; 



QL 



dm 
rf7 



- 2M 



(7 + 20 (to*)" - 13 (™^)^) + (5 + 8 (™^)' 



3 (68r - 141A/) + 48 (13r - 27Af) (rw"^)^ + 30 (14r - 29A/) (rM-^)" 



+ 48r {u^'f (3 + 5 (rw"^) 
3q^AP I r 



Ar^ + O (Ar^) 



(5.3c) 



11200rio V r - 2M 



1 + (ru"*)^ H (w'')^ 



X jlOr 

+ O (At 

3g^A/^ / 
^11200rio 1 



r-2M 



1 (^7 (™*) + 13 {ru^yj + (u'-y [1 + 8 {ru'f') 



At' 



(20r - 49M) - 8 (17r - 31M) (rw*)^ - 30 (14r - 29A/) (ra"^)"^ - 16r (m'')^ (^1 + 15 (ri/)' 



At" 
(5.3d) 



5r^ 



r - 2A/ 



5 + 28 {ru'l'y + 26 (ru"^)") + 4 (u'')^ (l + 4 (rit"^)^) 



At^ 



3u' 



(20r - 41A/) + 8 (17r - 35A/) {ru'^f + 10 (14r - 29A/) {ru'*')^ + 16r (u^'f (^1 + 5 {ru'^'Y^ 



At^ 



0(A7 



(5.3e) 



We see that the ^-component of the 4-acceleration is as would be expected for motion in the equatorial plane. We 
also find that the other three components, magL, maq^ and ma\^^ all have their leading order terms at O (At"*) and 
also have terms at order O (At^). 

We may also express this result in an alternative way by making use of the constants of motion [29j : 



r - 2M 



where 



= -Ve2_i_2Feff(r,e,0 



V,s{r,e,l)^-— + ^- — 



(5.4a) 

(5.4b) 
(5.4c) 



(5.5) 



is an effective potential for the motion. 

This gives a form for the equations of motion which will prove useful later (as a check on our result for Kerr 
spacetime): 



ma 



- 1 - 2Vcff (20 ~ 40 PMr^ + 50 l*r^ - 100 Af r"^ + 10 eV^ + 80 e^Pr'^) At'^ 



11200ri7 

+ (-60 fr^ + 255 fMr^ - 270 fM^r^ - 240 + 1008 l^Mr^ - 1056 1'^M^r^ - 180 l^r^ + 750 l^Mr 
- 780 fM'^ - 36 eV® + 81 e^Afr^ - 264 e^Pr^ + 552 eH'^Mr^ - 60 eH'r* + 90 eH^Mr^ + 16 e'r^ 
+240 e'^Pr^) Ar^ + O (At^)] (5.6a) 
ma" = O(AtS) (5.6b) 



11200r-i6 



[(20 - 40 Mr^ + 70 Pr'^ - 140 /^A/r^ + 50 Pr"^ - 100 PMr + 50 eV^ + 80 e^Z V^) At^ 



- x/e^ - 1 - 2V^s (60 - 135 Afr^ + 240 Pr^ - 528 ^^Afr^ + 180 Pr - 390 PM 
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+144 + 240 e^Pr^) Ar^ +0 (At^)] (5.6c) 



11200 (r - 2M) 

[(-10 + 20 Mr^ - 20 /V^ + 40 l^Mr^ + 50 ^^r^ - 100 l*Mr + 10 eV^ + 80 e^? V**) At'' 
- - 1 - 2Vcs (-36 + 81 Mr* - 120 Z^r''' + 264 pMr^ + 180 - 390 1'^M + 16 e^r^ 
+ 240 e^;^^) At^ +0 (At^)] (5.6d) 

c?m ^ ^^!^ [ (-5 r^^ + 10 A/r^ - 40 Pr* + 80 Z^Af - 50 l^r^ + 100 ^^Mr - 20 e^r*^ - 80 6^/2^4) ^ 4 
dr 11200ri4LV / 

+ Ve2 - 1 - 2yeff (12 - 27 Afr'^ + 120 Pr^ - 264 PMr^ + 180 ?'*r - 390 PM + 48 eV^ + 240 e^Pr^) Ar^ 

+0 (Ar«) ] (5.6e) 

These expressions give us the quasi-local contribution to the equations of motion for a particle following any geodesic 
path in a Schwarzschild background. It is also interesting to look at some specific particle paths. In particular, in 
order to check our results, we consider three cases for which results using other methods are already available in the 
literature [li,[l3^. 

1. A particle following a circular geodesic 

2. A particle under radial in- fall 

3. A particle under radial in- fall from rest 



1. Circular geodesic in Schwarzschild 

For a particle following a circular geodesic in Schwarzschild spacetime, we have the condition u'~ — along with 
the previous condition = 0. Furthermore, for a circular geodesic orbit u'^ is uniquely determined by the radius, r, 
angle. 6 and mass. A/ of the Schwarzschild black hole: 



> - 



1 / AI 



= — 



r\ r~ 3M 



(5.7) 



Substituting (|5.7p into (|5.2p . we also obtain the expression for m* in terms of the radius, r, and mass. A/ of the 
Schwarzschild black hole: 



3A/ 



(5.8) 



We can now use these values for the 4- velocity in Eqs. (|5.3ap - (|5.3ep . In doing so, the expressions simplify 

significantly and we obtain a result for the quasi-local contribution to the equations of motion for a scalar particle 

following a circular geodesic in terms of q, r and M alone: 

3fl2^/2(r-2A/)(20r3-81A//r2 + 54rAf2 + 53Af3^ 

ma^QL = — ^ 5 ^Ar^ + O (At*') (5.9a) 

^ 11200 (r-3A//)'rii ^ ' 

mo^L = 0{At^) (5.9b) 

J, 3q'^M^{r~2Mf{7r~8M) I a7 ^ 4 r. ^ 

maL = — ^ '—^ -\ — At'^ + 0(At6) (5.9c 

1120(r-3A./)'rio V r - 3A/ ^ ^ ^ ' 

, _ 3qHP{r-2M){7r-8M) I , 



ma*oL = — ^ -\ At^ + O At'' (5.9d) 

-At* + O (At^) (5.9e) 



1120 (r-3A//) r9 
dm 3q^AP (r - 2Af) (l3Af2 + 2Mr - 5r'^) 



dr 2240 (r - 3Af) 

Our results exactly match Eqs. (j6.4p - (|6.11bp containing our corrected version of the results of Anderson and 
Wiseman [l3|. 



^ Ref. [Ts^ also considers the case of a static particle, but in that ease the motion is non-geodesic, so we leave the analysis for later work. 
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2. Radial geodesic in Schwarzs child 



A particle following a radial geodesic in Schwarzschild spacetime has the conditions on the 4-velocity that u'^ = 
in addition to the previous condition — 0. We can therefore write Eq. 



r - 2M 



r - 2M 



{ury 



(5.10) 



This leaves us with one independent component of the 4-velocity, . Substituting for u*^ and into (|5.3ap - (|5.3ep . 
the expression simplifies significantly and we obtain a result which is only slightly more complicated than for the case 
of a particle following a circular geodesic. Our expression for the quasi-local contribution to the equations of motion 
is now in terms of the scalar charge, q radius, r, the mass, M and the radial 4-velocity of the particle, u^: 



maqL = 



mcQL = 



2 

- - 2M + r (u^'f^j IOt^u^'At* + (a9M - 20r + 16r (u^'f^j At^ + O {Ai 



'a*QL 

dm 
rf7 



3qHP 
11200r 
O (Ar^) 

O (At6) 

11200 (r - 2M)r9 



(5.11a) 
(5.11b) 
(5.11c) 



2M 



11200r 
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{u'-f lOr^u'^Ar^ + ( 20r - 49M - 16r (u^'f) Ar^ 4- O (At^) 'b.Ud) 



(5. lie) 



5r ( 5r - lOA/ + 4r {u''f ] Ar"* + u'' ( 60r - 123M + 48r {u^'f ] At^ + O (At 



From angular momentum considerations, we would expect a radially moving particle in a spherically symmetric 
background to experience no acceleration in the 9 or directions. This is confirmed by our results. 

3. Radial geodesic in Schwarzschild: Starting from rest 

In order to calculate the quasi-local contribution to the equations of motion for a particle under radial infall from 
rest, we begin with expressions (j5.11a[) - (|5.11ep for a particle under general radial infall. For the sake of simplicity, 
we will take r = to be the proper time at which the particle is released from rest at a point r = vq. We also assume 
that at proper time r = At, the particle will be at the radial point r, travelling radially inwards with 4-velocity u^. 
In this case, clearly expressions (|5.11ap - (|5.11ep hold. However, we must now obtain an expression for = {At) 
satisfying the initial condition (0) = 0. 

We begin with an expression for the radial 4-velocity in radial geodesic motion starting from rest at rp psj : 



(5.12) 



Assuming that the particle has not travelled too far (i.e. (r — ro)/ro "C 1), this equation may be approximately 
integrated with respect to r from tq up to r, giving an approximate expression for tq: 



ro =r 



MAt 



O (Ar^) = 



MAr2 

2r2 



O (At3) 



(5.13) 



Using this in (j5.12[) , we get an approximate expression for (we have taken the negative square root since we are 
considering infalling particles): 



A/A7 



O Ai 



(5.14) 



We can then use this expression in (jS.llap - (|5.11ep to give the quasi-local contribution to the equations of motion 
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for a particle falling radially inward from rest in Schwarzschild spacetime: 

iq^KP (r-2M\ 



mar 



j (20r - 39M) Ar^ + O (Ar*^) (5.15a) 



ma^L = O(Ar^) (5.15b) 

ma^L = 0(Ar'^) (5.15c) 

O (Ar*^) (5.15d) 

^ ' ^At^ + 0{At^) (5.15e) 



ma*QL 



dr 448r8 

As explained in Section IVIBl this result is in agreement with those of Ref . [I 



B. Kerr Space-time 

Next, we look at the much more complicated but physically more interesting Kerr spacetime. In Boyer-Lindquist 
coordinates, the line-element is 

,2 A 2Mr\ ,2 4aMrsin2 6' , S,. v. ,.2 2Mr{r^+a^)\ . , . , 
ds^ = - il — j dt^ dtdc/} + —dP + Sd6|2 + f A + ^— '- j sin^ Odcj)^ (5.16) 

where 

S = P + a^cos^e (5.17) 
A = - 2Mr + a^ (5.18) 

Although it is possible to use a computer algebra package to obtain a totally general solution for the self-force in Kerr 
spacetime, the length of the resulting expressions are too unwieldy to be of much benefit in printed form. For this 
reason, we will only focus here on equatorial motion, so that 9 = keeping in mind that more general results may 
be obtained using the same methodology. 

As a result of the axial symmetry of the metric, motion in the equatorial plane in Kerr spacetime will always 
remain in the equatorial plane, i.e. = 0. Additionally, the motion is parametrized by two quantities analogous to 
the Schwarzschild case: the conserved energy per unit mass, e and the angular momentum per unit mass along the 
symmetry axis, I. We can therefore write the three non-zero components of the 4- velocity in terms of these conserved 
quantities [29j : 



u 



u 



dr 
d^ 

dr 
dt 
d7 



where 



v/e2-l-21/eff(r,e,0 

r - 2M\ , 2 A/a ' 
1 + e 



2 I 2 
r + a 



2Ma' 



2Ma, 
e 1 



Vce{r,e,l) = h 



M P - a2 (e2 - 1) M {I - aef 



2P 



(5.19a) 
(5.19b) 

(5.19c) 
(5.20) 



is the effective potential for equatorial motion. 

Substituting these values into Eq. (|4.ip and computing the relevant components of the Hadamard/DeWitt coef- 
ficients, we arrive at a lengthy but general set of expressions for the equations of motion for a particle following a 
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geodesic in the equatorial plane in Kerr spacetime: 



11200ri4 



1 - 21/eff 



(20 Pr^ - 40 l^Mr^ + 50 l^r"^ - 100 1'^Mr^ + 10 e^r^ + 80 e^Z V) 



3^2 



400 e^^Mr^ - IGOe^^r") 



(150 r 



2^2 



FSOeZA/r^-eOOeZ-r 
500 + 100 e^r^ - 40 e^Mr^ + 1500 e^^V^ - 600 e^PMr^ + 80 e^r^) 



(-300 eZr^ - 2000 - 1400 e^lr^ + 400 e^lMr^) 
(150 eV^ + 3000 e^l'^ + 450 eV - 100 e^'Mr^) 



-a^ (2000 e^Z) + (500 
3g2Af2 



11200ri8 



(60 Pr'^ - 255 PMr^ + 270 Z^AfV'^ + 240 Pr^ - 1008 ^''Mr^ + 1056 PM'^r^ + 180 Pr^ 



750 ^^A'/r^ + 780 PM^r + 36 eV^ - 81 e^Mr 
90 e^Z^'Mr'^ - 16 eV - 240 eH'^r'^) 

- 432 e^r^ + 1176 elMr^ ~ 540 e/M^r^ - 2760 e/V + 7728 ePMr^ 



552 e^Z^Mr^ + 60 e^^V 



4224 eZ^'M^r^' - 2160 e/V^ 



+ 6600 ePMr'^ - 4680 eZ^Af V - 368 e^Zr^ + 1104 e^/Af + 1680 e^/^r^ + 360 e^PMr^ + 480 e^/r^) 



-(660? 



2^5 



-921e^A/r6 



270e2AfV^ + 6720e2/V- 



21750e^^'^A/r2 



1368 PMr'^ + 3000 ? V - 6120 Z'^Af + 2100 Pr - 4200 Z^Af + 456 eV^ 

2/^^ - 17136 eH^Mr"^ + 6336 e^PM^r^ + 3120 e^/V^ 
104 eV - 552e'^A/r^ - 5400 6^^?^^ - 5A0e^fMr^ - 240 eV) 
13600 V + 24480 ePMr^ - 8400 ePr + 25200 ePM ~ 6120 e^lr^ 



U700 e^pM^r^ 



4224 eHM'^r^ + 5520 e^Z + 36000 eH^Mr'^ 



15600 e^Z^AfV 



5520 e^Zr^ 



( - 1800 eZr^ + 2736 elMr^ 
+ 15120 e^ZAfr'' 
+ 360e5;Afr'') 

a''(720 Pr^ + 2800 Pr + 1140 eV^ - 1368 e^Mr^ + 22800 e^ZV^' - 36720 eH'^Mr'^ + 10500 e^^r 



1056e''A/V3 - 16380 e^'ZV^ - 32250 e'^/^Afr^ 

'( - 1440 elr^ - 11200 el\ - 16800 eHr^ + 24480 e^/Afr^ + 84000 e'PM + 13440 e^/r^ 
+ 15000 e^ZA/r^ - 4680 e^/AfV) 

' (720 eV^ + 16800 e^/ V + 4600 - 6120 e^Mr'^ - 10500 e''/ V - 63000 eH'^M - 3720 e V 
- 2850 e^Afr^ + 780 e^A/V) 



63000 e^Z'' AT + 1920 e^r"^ - AIM e*Mr* 
11700 eH'^M^r - 1860 e^r^ - 90 e^Mr*) 



-11200 e^Zr + 8400 e^lr + 25200 e'^ZAf) 



(2800 eV - 2100 e'^r - 4200 e^Af) 



Ar^ + O {At^) 



O At" 



(5.21a) 
(5.21b) 



11200 (r2 -2A/r + a2)ri6 
+ a(60 er' 



( - 20 Ir^ + 80 /A/r^ - 80 IM'^r^ - 70 Pr" + 280 /^Afr^ - 280 PM'^r^ 



200 PM^r^ 



hOeHr 

2 I 

r 
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lOOe^ZA/r^ - 80 6^/^ + IQO e^^Mr^) 



160 eAfr^ + 80 eM\^ + 570 el\" - 1560 e^^A/r^ + 840 ePM^-'^ 



600 e;^-'^ 



1700e/''Afr^ 
- 170;r^ + 340/A/r5 



1000 eZ^AT V + 40 e V 



100 at/ 



700 ZV^ + 1400 Z^Afr^ 



160e^;V -480e^;2Afr'^) 
500 Z'V^ + 1000 PMr - 890 e^^r^ 



2280 eHMr^ 



MOe^lM'^r^ 



1500 e^^V^ + Ad,00 e^PMr^ - 2000 e^P M'^r'^ - 80e^Zr^ + 480e^;Afr^) 



(210 er" - 340 eMr^ + 2250 e/-r 



2000 e^V 



1000 e-'Mr" + 280 e^AfV'^ 



- 4200 eZ^A/r^ 
1400 e^ZV^ - Q200 e^fMr^ 



5000 e^^'A/r + 390 eV 
2000 e-'^Z^A/V^ - 160e5A/r^) 
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a^{-150lr^ - 500 l^r^ - 2400 e^lr"^ + 4200 e^lMr^ - 3000 e^l^r^ + 10000 e^/^Mr - 450 eHr'^ 
+ 3800 eHMr^ ~ 1000 eHM^r^) 

(150 er^ + 1500 efr^ + 850 eV'^ - 1400 e^Mr^ + 2000 e^fr'^ - 10000 e^P A/r - 900 e^Mr^ 
+ 200e^AI^r^) 

( - 1500 e^lr^ - 500 e'^Zr^ + 5000 eHMr) 



a^(500eV - 1000 e^Mr 
3qHP 



v/e2 - 1 - 21/e 



cff 



I - 255 IMr^ + 270 ^A/V^ + 240 1 



3^5 



11200 (r2 -2Afr + a2)ri6 

- 1008 /^Afr'' + 1056 PM^r^ + 180 - 750 l^Mr^ + 780 l^M^r + 144 e^/r^ - 288 e^lMr^ 
+ 2'iOe^Pr^ - 480 e^l^Mr"^) 

+ a{~ 192 er^ + 528 eMr^ ~ 270 eA/V^ - 2040 el'^r^ + 5688 efMr^ - 3168 efM^r^ - 2160 e/V 
+ 6240 eZ^/r^ - 3900 el'^M^r - 128 eV^ + 288 e^A/r*' - 480 e^Pr^ + 1440 e^PMr^) 

+ (660 ?r'5 - 1368 IMr^ + 3000 ? - 6120 PMr^ + 2100 Pr - 4200 Z^A/ + 3360 e^lr^ - 8352 e^/Afr'^ 
+ 3168 e^lM^r^ + 5400 e^Pr^ - 17460 e^PMr^ + 7800 e^PM^r + 240 e-^/r^ - 1440 eHUr^) 

+ a^{- 840 er^ + 1368 eAf r'' - 10080 el\^ + 18360 e/^A/r^ - 8400 ePr + 21000 ePM - 1560 eV^ 
+ 3672 e^Mr^ - 1056 e^A/V^ - 5040 e^Pr^ + 22440 e^PMr^ - 7800 e^fM\ + 480 e^A/r'') 

+ 0^^(720 + 2800 + 11160 e^lr^ - 18360 e^/Afr^ + 12600 e^Pr - 42000 e^/^A/ + 1620 eHr^ 

- 13710 eHMr^ + 3900 eHM^r) 

+ a^{- 720 er^ - 8400 efr - 4080 eV^ + 6120 e^Mr^ - 8400 e^fr + 42000 e^?^ Af + 3240 e^Mr"^ 

- 780 e'^A/^r) 



a 



e^lr + 2100 e^'/r - 21000 eHM) 



ma* = — 



(- 2800 eV + 4200 AT 
3qH'P 



At^ + O (Ar^) (5.21c) 
(10 er^o - 20 eMr^ + 20 ePr^ ~ 40 e/^Afr^ - 50 e/ V + 100 e^^Mr'^ 



11200 (r2 -2A/r + a2)ri6 

- lOeV^o - SOeH^r^) 

+ a{- 60 Ir^ + 160 /A/r^ - 80 IM^r^ ~ 150 l^r^ + 440 Z^Af - 280 /^A/^r"* + 100 /^A/r^ 

- 200 l^M'^r^ + 180 eHUr'^ + 600 e^ZV^ - 240 e^PMr^ + 160 e^Zr®) 

+ a2(llOer^ - 240 eA/r^ + SOeAfV + 620 eZV^ - 2040e/2A/r^ + MOel'^M'^r* - 550 eZV 

- 1200eZ''Mr3 + 1000 eZ'' A/ V - 30 eV - UOe^AV - 1580 e^ZV^ + 120 e^fMr^ - 80 eV®) 



a 



3 210 Ir^ + 340 IMr^ - 650 + 1400 PMr^ + 1000 Z^A/r - 750 eHr^ + 2760 eHMr^ 



- MOeHu'^r^ + 2600 e^^V^ + 3800e2;^A/r3 - 2000 eH^ M"^ r"^ + 1560 e^^r^ + 80e''/Afr^) 
a^(250er'^ - 340eArr^ + 2100 e/V - 4200 efUr^ ~ 500 e/V - 5000e/''Afr + 280 eV - 1160e^Afr^ 

+ 280 e^M'^r^ - 4500 e^^V - h200 eH'^Mr^ + 2000e^Z2A/V - 530 - 60e^Afr5) 
( - 150 - 500 - 2250 eHr'^ + 4200 eHUr^ + 2000 e^ZV^ + 10000 eH^Mr + 3400 eHr'^ 
+ 3300 eHMr^ - 1000 eHu'^r'^) 

(150 er"* + 1500 efr'^ + 800 - 1400 e^Mr^ - 3000 e^^V^ - 10000 e^P ATr - 950 eV^ 

- 800 e^ATr^ + 200 e^M^r^) 

( - 1500 e^/r^ + 2000 eHr^ + 5000 e''/ Afr) 
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*(500eV - 500 eV - 1000 e^Mr) 



At' 



11200 {r^ -2Mr + a^): 



.16 



1 - 2v; 



off 



36 er^ + 81 eMr^ - 120 elV + 264 el'^Mr^ 



dm 
d7 



+ 180 V - 390 e^^'A/r'' + 16 + 240 e^/ V^) 
+ (240 /r^ - 648 IMr^ + 270 ^MV^ + 720 l^r^ - 2040 l^Mr' + 1056 /^MV^ - 360 l^Mr'^ + 780 Z^A/V 

+ 240e2/r^ - 816 e^/A/r^ - 2160 e^Pr'' + 1080 e^PMr'^ ~ 480 eHr'^) a 
+ ( - 456 er^ + 921 eMr^ - 270 eA/V^ - 3360 el\^ + 8784 e^^Af r'' - 3168 efM^r^ + 2280 V 

+ 4290 el'Mr^ - 3900e/'^A/V - 104 eV + 552e3A/r^ + 5640 e^/V^ - 900 e^f Mr' + 240e'V^)a^ 
+ (960 Ir^ ~ 1368 IMr' + 3520 - 6120 Z^A/r^ - 4200 l^M + 4560 e^lr^ - 11448 e^/Af 

+ 3168 e^lM^r^ - 10560 e^ZV^ - 13560 e^/^A/r^ + 1800eH^M'^r ~ hh20eHr' + 120 e^ZAfr'')a^ 
+ ( - 1140 er^ + 1368 eMr' - 11640 e^V^ + 18360 el^Mr^ + 2100 e^V + 21000 el'U - 1920 eV^ 

+ 4704 e^Afr-* - 1056 e'M'^r' + 18000 e^/^^s + 18540 eH^Mr^ - 7800 e^l^M^r + 1860 e^r^ 

+ 90e5Afr^)a^ 

+ (720 Ir^ + 2800 l^r + 12720 e^Zr^ - 18360 eHMr^ - 8400 e^/^r - 42000 e^l^M - 13440 eHr^ 

- 11760 e'lMr^ + 3900 eHM^r)a'^ 
+ ( - 720 er^ - 8400 e^V - 4600 e^r^ + 6120 e^Mr^ + 12600 e^ZV + 42000 e^PM + 3720 eV^ 

+ 2850e^Afr2 - 780 e^A//2r)a^ 
+ (8400 e^lr - 8400 eHr - 21000 eHM)a'' 

+ ( - 2800 e\ + 2100 + 4200 e^Af) a® Ar^ + O ( Ar^ 

3q2Af2 



(5.21d) 



11200ri5 



(-5 + 10 A/r^ - 40 fr"^ + 80 fMr' - 50 ^''r^ + 100 1'Mr^ ~ 20 e^r^ - 80 e^fr^) 



. (240 elr^ - 160 elMr' + 600 el^r^ - 400 ePMr^ + 160 eHr^^ 

2 (-20 - 300 - 500 1'r - 200 eV^ + 80 e^AIr' - 1500 e^fr^ + 600 e^fMr^ ~ 80 e 
^ (600 elr^ + 2000 el\ + 1400 eHr^ - 400 e^Mr^) 
' (-300 eV - 3000e2;2r - 450 eV + lOOe^A/r^) 



+ (2000 eV/) - (500 e'r) 
3q^M^ 



At" 



11200ri5 



- 1 - 214ff 



12 - 27Afr^ + 120 - 264 ^^Mr^ + 180 - 390 1'Mr 



+ 48eV + 240e2;V) 



( - 720 elr' + 528 elMr^ - 2160 e/ V + 1560 el^Mr - 480 e^/r"*) 

- 1080 + 2100 + 600 e^r' - 264 e^Mr^ + 5400 e^/V^ - 2340 e^l^Mr + 240 eV) 
^( - 2160 eZr^ - 8400 eZ^ - 5040 e^Zr^ + 1560e^/Afr) 
' (1080 eV^ + 12600 e^P + 1620 e V - 390 e'Mr) 



a^(60 



0^(8400 6^?) +a6(2100e'' 



At'5 + O (At*') 



(5.21e) 



It is straightforward to see that these results exactly match Eqs. (|5.6a[) - (|5.6c|) in the limit a ^ 0. 
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1. Release from instantaneous rest 

We now consider the motion of a particle which is released from rest relative to an observer at spatial infinity. This 
is analogous to the case of radial infall from rest in Schwarzschild (Section IV A 3|) and we proceed with the calculation 
in exactly the same way. Imposing the initial conditions u^{r = tq) = and u'^{r = tq) = 0, Eqs. (|5.19ap and (j5.19b[) 
can be solved for the constants of motion: 



ro - 2M 



I 



2Ma 



ro 



ro V ^0 - 2M 

We then substitute these into Eq. (|5.19ap and integrate to get an approximate expression for rg in terms of r 



(5.22a) 
(5.22b) 



ro = r + (Mrl - (P - (e^ - l)) ro + 3M (/ - aef) ^ 
V / Zr^ 

This, along with the equations for e and I then allow us to express the 4-velocity as: 

V2 _ 2j.M + a^' 



u = — 



r3 (r - 2M) 



MAt 







= - 



r - 2M \ 2Ma 



r 



2Ma r~2M 



1 

A 



r I r y r — 2AI r 

2Ma^ \ /r-2M 4A'Pa^ 



r 



2M 



(5.23) 

(5.24a) 
(5.24b) 
(5.24c) 

(5.24d) 



We then substitute these values along with the components of the Hadamard/DcWitt coefficients into Eq. (|4.1 
and project orthogonal and tangent to the 4-velocity to get the equations of motion: 



ma = — 



11200 (r2 - 2 Mr + a^f (2 M - rf r^^ 



- 4a^ (252 - 160 rM^ - 162 r^M + 95 r^) (2 M - r)° r 



- (39M-20r) {2 M - rf r'^^ 



6 „12 



12a^ (848 Af^ - 640 rA/^ - 1404 r^Af^ + 960 r^Af^ + 291 r^A/ - 185 r^) (2 M ~ rf r^ 



Aa" (12288 A/' - 10240 rA/° - 42720 r^A/^ + 32640 r^Af* + 20844 r^'M^ - 14880 r^'A/^ - 2292 r'' AT 



„4 njS 



y.5 i\,r2 



4 ^6 



h 1545 r^) (2 M - r) 

(98304 A/^ - 81920 rA/^ - 811008r2Af^ + 655360 r^A/^ + 758592 r^A/^ - 583680 r^M'' 
- 187920 r^ A/3 + 138240 r^A-/^ + 13347 r^M - 9420 r^) (2 A/ - r) 



3^3 



12a^" (8 A/2 - r^) (15872 A/^ - 12800 rA/^ - 30272 r^M^ + 24000 r^A/^ + 10704 r^M^ - 8200 r'^A/^ 
- 921 r^M + 675 r^) (2 A/ - r)^ r^ 



lOa'-' (2 A/ - r) (7296 A/^ - 5760 r A/* - 4704 r^A/^ + 3680 r^A/^ + 489 r^'M - 370 r") (8 A/ 



ma^ 



+ 100ai'^ (9A/-7r) {8 - r^f 
O (Ar6) 

Sq^M^a 



112 (-2 A/r + r2 + a2)3 (2 A/ - rf^^ r^V^ 
+ a^r^ (16 A/2 - 7 r^) (2 A/ - rf 



Ar^ + O (Ar^) 



r^ (2 A/ - r)^ 



(5.25a) 
(5.25b) 
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+ 16 a^r^ {4:M^~6 M^r^ + r^) (2 M - rf 
- 5 a^r (2 M - r) (-3 + 8 A/^) (8 A/^ - r^) 



5 a* 



ma 



(5.25c) 



56 (r2 - 2 A/r + a2)^ (2 M - r)^/^ r43/2 

„2 



(2A/-r)^ 



dm 
d7 



(40 Af^ - 16 rA/2 _ 30 r^M + 9 r^) (2 M - r)^ 



a^e 



(256 Af^ - 64 rA/^ - 560 r^A/^ + 128 r^A/^ + 138 r^M - 31 r^) (2 A/ - r)* 

- 2a^r3 (256 Af^ - 1664 AfV + 224 A/V + 1100 A/V - 184 AT V - 147 Mr^ + 27 r^) (2 A^ - rf 
+ c?r^ (6144 A'f - 10752 A-f + 1024 A/V + 3680 MV^ - 496 A//V^ - 324 A/r^ + 51 r^) (2 M - r) 

- 5 a^" r (2 M - r) (-5 + 36 r^M + 24 r^Af ^ _ 272 r^Af ^ + 384 M'^) (8 Af^ - r^) 

V 2 



5ai2 (^32 m3 _ g ^2^^ 



3(72 m2 



Af 2 _ r' 



At^ + O (Ar^) 
r^^ (r-2Af)^ 



(5.25d) 



448r20 (r2 - 2 A/r + a2)'' (2 Af - rf 
+ 16 aV^^ (r^ - 2 A/^) (r - 2 Af )^ 
+ 12 aV^ (7 - 40 r^A/^ + 32 M^) (r - 2 Af )^ 
+ 8 a^r^ (27 r*^ - 276 Af V + 672 Af ^r^ - 256 M*^) (r - 2 Af)^ 
+ aV^ (r^ - 8 A/2) (-512 Af ^ + 3264 Af^r^ - 2232 AfV + 309 r^) (r - 2 Af)^ 



12 a^^r^ (21 ~ 96 r^A/^ + 64 M^) (r - 2 A/)^ (r^ 



8A/^ 



2\3 



+ 10 a^^r (r - 2 Af ) (-24 Af ^ + 11 r^) (r^ - 8 A/^) 



+ 20 a" {f -~8M^y 



0(A7 



(5.25e) 



Again, it is straightforward to see that these results exactly match Eqs. (|5.15ap - (|5.15cp in the limit a — + 0. 



VI. COORDINATE CALCULATION 

Using a Hadamard-WKB expansion, Anderson and Hu [3l|^ calculated the tail part of the retarded Green's function 
in Schwarzschild spacetime in terms of the Schwarzschild coordinates up to O ^(a: — x')^^ in the separation of the 
points. They then used their results as a basis for the calculation of the entire self- force on a particle held at rest 
until a time t = 0, at which point it is released and allowed to fall freely [l^. Due to the spherical symmetry of the 
Schwarzschild spacetinie, the resulting motion is radially inwards. 

The results of Ref. [3l| were also used by Anderson and Wiseman [13] to calculate the self-force on a particle in 
Schwarzschild spacetime under two other circumstances: a static particle and a particle following a circular geodesic. 
Unfortunately, although the technique they describe is valid, there are some errors in the results. For this reason, 
and in order to extract the radially infalling part from the results of Ref. [iBl , we reproduce those calculations here 
with corrections. We also go one step further and calculate the equations of motion from the self- force. Although the 
static particle result also contains errors, we will not consider it here as it corresponds to non-geodesic motion, which 
we have considered, but leave for later work to avoid confusion. 

We follow the method prescribed by Anderson and Wiseman (l7j and begin with the same expression we had for 



It is important to use the corrected expressions for the n^jfe as given in the erratum to this paper. 
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the covariant calculation: 

/ql = -9'/ V-v{x,x')dT' (6.1) 

Next, instead of using a covariant expansion for v (x, x') as we did previously, we now use the coordinate expansion 
given in Ref. j3ll|'^: 

oo 

v{x,x')= v,,k{t-t'ficos{-f)-iy{r~rf (6.2) 

The coordinate 7 is the angle on the 2-sphere between x and x' . For the spherically symmetric Schwarzschild 
spacetime, we can set ^ = f without loss of generality, so that 7 = 0. We then take the partial derivative with 
respect to the coordinates r, 0, and t and convert from those coordinates to proper time r. Next, we integrate over r 
as we did previously in the covariant calculation and project orthogonal to and along the 4-velocity to obtain a final 
expression for the quasi-local contribution to the equations of motion. This coordinate transformation is dependent 
on the particle's path so it is in this way that the motion of the particle affects the self-force and equations of motion. 
We now work through this calculation for two of the particle paths investigated previously in the covariant calculation. 



A. Circular Geodesic 



For a circular geodesic, the geodesic equations can be integrated to give 

r-r' = 0, e~e' = 0, cj)-(b'^-J ^I^A t-t'), t-t' 



r-3M 



(r-r') 



(6.3) 



We now substitute Eq. (|6.2p into (|6.ip . take the partial derivative, then use the above to express the coordinate 
separations in terms of proper time separations. Finally, we do the easy integration over r' and project orthogonal 
and parallel to the 4-velocity to obtain an expression for the quasi-local contribution to the equations of motion for a 
particle following a circular geodesic in Schwarzschild spacetime. The radial component of the quasi-local mass times 
4-acceleration is 



where 



[5] = - 



[7] = - 



ma£jL = ma^QL [5] At' + ma'^^ [7] Ar^ + O (Ar^) 



Sq^R'P (r - 2M) {20r^ - 8lMr^ + 5AAPr + 
11200 (r-3M)%ii 
q^M^ (r - 2M) 
188160 (r - 3M)^ 

(560r^ - 6020r^M + 21465r^M2 - 26084r2M^ - 5281rM^ + 21510M^) 



The 6 component of the quasi-local mass times 4-acceleration as would be expected. 
The (p component of the quasi-local mass times 4-acceleration is: 



ma, 



QL 



ma, 



QL 



[4] At^ 



[6] Ar^ + O (Ar« 



where 



ma^L [4] 



3q^AP {r - 2M) {7r - 8M) 
1120 {r-3Mfr^o 
q^AP [r ~ 2M) 



M 



r-3M 



13440ri3 (r - 3M) \l r - 3M 



M 



(6.4) 
(6.5a) 

(6.5b) 

(6.6) 
(6.7a) 



- (I26r* - 1129r3A/ + 3U7M\^ - A19MI^r + 1623M^) (6.7b) 



* There are several different conventions in the literature for the Hadamard form of the retarded Green's function. For the present work, 
we use a convention consistent with Ref. [23l throughout. However, the convention of Ref. |31| differs by a factor of ^ in the coefficient 
of V {x, x'). Fortunately, our convention for the scalar charge, q also differs by a factor of V47r, with the result that our covariant V (x, x') 
is exactly equal to the coordinate v{x,x') of Anderson and Hu [31II . Furthermore, Ref. |31| originally indicated a factor of in the 
coefficient of v (x,x'). It has subsequently been confirmed with the authors that this should be -tj:. 
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The t component of the quasi-local mass times 4-acceleration is: 



ma^T. = ma*nT, \4] At"^ + ma*QL [6] Ar^ + O (Ar^) (6.8) 



where 



ma. 



1120(r - 3A/) r9 \l r - 3M 



41 ^ (6^9a) 



ma. 



Ql[6] = _ , — ^_^3 '^^^ (l26r^-1129r^M + 3447r^M^-4193rA/^ + 1623M^) (6.9b) 



13440 (r - 3A/)-^ V r - 3M 
Finally, the quasi-local contribution to the mass change is 



where 



dm dm . ^. ^ . dm . c „ / . on , „n 

— = — 4 At^ + — 6 Ar^ + O (Ar^) (6.10) 
dr dr dr 

dm , , 392 A/2 _ 2M) (UAP + 2Mr - Sr^) 

— 4 = — ^ '-^ 5 ^ 6.11a 

dr 2240 (r - 3Mf 

dm , , q^M^ {r - 2M) (651 - 344 M^r - 261 M^r^ + 189 Mr^ - 28 r^) 

— — 6 = H (6.11b) 

dr 6720 (r - 3Af )^ 

Direct comparison between these results for the equations of motion and those of Ref. [l3l for the self-force are not 
immediately possible. However, it is straightforward to project the results of Ref. 17[ orthogonal to and tangent to 
the 4- velocity. In doing so, we find that both sets of results agree in their general structure, but differ in the numerical 
coefficients involved. We suspect a minor algebraic slip may have led to the errors in the results of Ref. [T7| . 



B. Radial Geodesic: Infall from rest 

The case of a particle following a radial geodesic has been investigated by Anderson, Eftekharzadeh and Hu [l6[ . 
They consider a particle held at rest until a time i = 0. The particle is then allowed to freely fall radially inward. In 
such a case they have shown that, provided the particle falls for a sufficiently short amount of time, it is possible to 
obtain an expression for the entire self-force, rather than just the quasi-local part. Unfortunately, direct comparison 
with our result is not immediately possible as their results give the full self-force, rather than just the quasi-local 
contribution from the radial infall phase. However, as it is straightforward to calculate quasi-local contribution using 
the method they prescribe, we have done so and found agreement with our results once the corrections to Ref. |3l| 
mentioned earlier in this section are taken into account. 

As we are only interested in the leading order behavior for comparison to our own results, we take a slightly more 
straightforward approach to the calculation. From the arguments in Section IV A 3[ it is clear that 



--'^'--$(-'--")' 0-e' = Q, 0-0' = O, t-t'^^-^^{r^r') (6.12) 

As in the circular geodesic case, we now substitute Eq. (|6.2p into (|6.1|) . take the partial derivative, then use Eq. 
(|6.12p to express the coordinate separations in terms of proper time separations. Finally, we do the integration 
and project orthogonal and parallel to the 4-velocity to obtain an expression for the quasi-local contribution to the 
equations of motion for a particle starting at rest at r = rg , t = and subsequently falling radially inwards for a 
proper time At to a radius r 

r - 2M 







3q^M^ / 
11200?'!" V 


ma^L 


= 


+ (At^) 


ma^L 


= 


+ (At^) 




= 


+ (At6) 


dm 




3g2M2 /r 


d7 




448r8 V 



(20r - 39M) Ar^ + O (At^) (6.13a) 

(6.13b) 
(6.13c) 
(6.13d) 

^^^^ Ar^ + O (Ar^) (6.13e) 



This is in exact agreement with our covariant results in Eqs. (jS.lSap - (jS.lSep . 
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Local Truncation Error 



As our results differ from those of Ref. [17[, it is desirable to also update the plots of the relative error given there. 
We estimate the local fractional truncation error, as Anderson and Wiseman have done, by the ratio between the 
highest order term in the expansion ( O (At"), say) and the sum of all the terms up to that order. 



/ql ["] 



(6.14) 



For the circular geodesic motion investigated in the present work, we plot in Fig. [l]the fractional local truncation 
error for the equations of motion as a function of Ar for r ~ 6Af , r ~ lOM, r — 20M and r = lOOAf, as was done 
previously for the self- force in Ref. [l7| . 





At(M) 




At (M) 



r=6M r=10M 



-r=20M 



-r=100M 



FIG. 1: Local truncation error in the equations of motion for a particle following a circular geodesic in Schwarzschild. 



VII. CONCLUSION 



The self-force is determined by the entire past history of the particle, while the quasi-local contribution is the 
portion of the self-force resulting from the particle's recent past. In this paper, we have calculated the quasi-local 
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contribution to the self-force and to the equations of motion for a scalar charge moving in a curved background 
spacetime. We have determined these quantities up to fifth order in At, the proper time distance backwards along 
the particle's path to which the quasi-local contribution applies. In addition to this result, which applies to a general 
spacetime, we also present the same result for vacuum spacetimes and show how the expressions simplify significantly. 

We then looked at some specific examples of spacetimes and particle motions. We looked at the Schwarzschild and 
Kerr spacetimes and at various geodesic motions. This allowed us to compare our results to existing work and to 
verify the validity of our approach. It also allowed us to find some errors in some of the existing work. 

While Anderson, Eftekharzadeh and Hu have successfully calculated the quasi-local contribution to the self-force 
to a much higher order (in At) than we have done here, their results apply only to a specific choice of spacetime 
and to a particular particle motion. On the other hand, the covariant method we describe only provides terms at the 
leading two orders for massless fields in vacuum spacetimes, but does so for a general vacuum spacetime and general 
(geodesic) particle motion. In this sense, the two methods are seen as complimentary, each proving as a useful check 
on the others results. 

We have chosen to restrict the present work to geodesic motion. This was desirable primarily in order to avoid 
confusion. Non-geodesic motion adds some extra complexity to the covariant calculation due to there being two 
different proper times: the particle proper time and the geodesic proper time. We have made significant progress in 
dealing with this issue and intend to present our solution to the problem at a later date. 

Obviously, the scalar case as presented here is of less immediate physical interest than the gravitational case. The 
biggest obstacle to achieving results in the gravitational case will be to calculate the relevant DeWitt/Hadamard 
coefficients. While this has previously been done for vacuum spacetimes up to 0(cr'^/^) by Anderson, Flanagan and 
Ottewill [ll], it is desirable to take the calculation to higher order to allow the matching point At to lie further 
into the particle's past. It would also be desirable, but not necessary, to obtain expressions for a general spacetime, 
rather than restricting ourselves to only vacuum spacetimes. Achieving these goals would be possible using a method 
similar to that of Folacci and Decanini [2^ . for example. However, it is probable that the expressions would be even 
longer for spin-2 (gravitational) than for spin-0 (scalar) fields. Some hope lies in the fact that we are most interested 
calculating V{x,x') for massless fields in 4-dimensional vacuum spacetimes. There is also hope that work by Fulling, 
et al. [1^ and more recently by Decanini and Folacci [13] to establish a linearly independent basis for the Riemann 
polynomials will simplify the calculation somewhat. Once V{x,x') has been calculated to sufficiently high order, it 
will be straightforward to apply a method analogous the the one presented here to calculate the self-force. 
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APPENDIX A: CALCULATING COVARIANT HAD AMARD/DE- WITT COEFFICIENTS TO S"' ORDER 

IN THE SEPARATION OF THE POINTS 

Folacci and Decanini [22| have calculated the DeWitt and Hadamard coefficients of the Feynman propagator in a 
general spacetime up to O (tr^) in the separation of the points (i.e. fourth order in the geodesic separation of the 
points). In particular, we are interested here in the Hadamard coefficients for the function V {x,x'), which appears 
in the Hadamard form of the Feynman propagator. 

This may be extended to O (cr^/^) relatively easily by taking into account the symmetry V{x,x') = V{x',x) of 
the Green's function. Repeatedly taking symmetrized covariant derivatives of both sides of this equation and taking 
coincidence limits, we arrive at an expressions for the odd order terms in terms of all the lower (even) order terms 
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(Al) 



Recursively applying this identity to eliminate all odd order terms, we get expressions for the coefficients of the 
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terms of order cr 2 : 

15 5 , 

Vabcde = -^V.(^abcde) " -^'"{ab-cde) + ^'^ iabcd;e) (A2a) 

1 3 

Vabc = -jV-(abc) + 2"(a'';c) (^^b) 
Va = ^W;a (A2c) 

These expressions are equally valid for any symmetric bi-tensor, so we may use it for both the calculation of higher 
order Vnai...ap as used in (|3.1[) and higher order Vai...ap as used in (|3.7p . 

APPENDIX B: VACUUM EXPRESSIONS FOR THE HADAMARD/DEWITT COEFFICIENTS 

The expressions given in Ref. [23| are unnecessarily long and unwieldy for the purposes of the present work. We 
are currently only interested in massless fields in vacuum spacetimes (mfidd = 0, Raf3 = 0), so the majority of the 
terms vanish and we are left with much more manageable expressions for the w,iai...ap- They are: 

vo = (Bla) 
Voa ^ (Bib) 

VOab = ^^9abl (Blc) 
"VQabc = -^g{abl;c) (Bid) 



9 9 1 

525 (aklb^'-' ^Qg"-' (o^\p<yr\b:cd) ('^I'^l'' <Ap\d)\r 



1 9 2 

(aA^\P'^T\c:d) - Y575 L.p(a|^|bO|^|c|,,|d) - -^L- (a'-^lprl b^-^klclKld) 

8 - 4 

^Y575'"'"'V'"|p| lr|b'^k|c|K|d) - Y^C'^'^^C'i^^i bCi^idKid) (Ble) 

vi = —I (Blf) 
720 ^ ' 

= (Big) 
1 . 1 29 . 

^^''^ ^ ^1400*"'"^^'''^'"^^^°'''^'' ^ 1575^'"^^°'"'^'"^^'' ^ 25200^'"^^'''^'"^"'''''"''^ 



1680 po-rb 1344 ;a^PCTTK;b + ^^g'-- '-^paa'-^rK 



Ah 



•-^ ^para'^KX b + 1 onnn '-^ptrra'-^KA fc K^'^"') 



1800 ''^ " 18900 

^ /" ' I — \r^paTK ^ 

6720 8960 



"720 8960 



9072 " 18440 

where Cobcd is the Weyl tensor and / = CafSjgC"^'^^ . 

Moreover, we have chosen to work with an expression for V{x, x') as given in p.7p . rather than the expression used 
in Ref. [1^ and Eq. (|3.ip . The two sets of cocfhcients may be related by using 

2a = a'^a.c. (B2) 
to rewrite p.ip in terms of a series expansion in powers of ct'". We can then equate, power by power, the coefficients 
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V 


— 7Jn 

— ty(j 




Va 


= VOa 




Vab 


= VOab + VlQab 




Vabc 






Vabcd 


= VOabcd + ^VK^abdcd) 


+ Gv2gi^ab9cd) 




1 5 


5 


'^abcde 







of powers of tr'" to the corresponding terms in (|3.7p to give: 

(B3a) 
(B3b) 
(B3c) 
(B3d) 
(B3e) 

(B3f) 

Filling in the values for the Vnai...ap we see that a large number of terms either cancel or combine to give: 

V = (B4a) 

Va = (B4b) 

Vab = (B4c) 

Vabc = (B4d) 

1 2 1 

Uabcd — 280 (° b\;a\ \p\c\c7\d) 315 '^pia\T\b'^\a\c\K\d) + IQ5 ^ (a h"-' \p\c'-'\TK.a\d) 

+ -^C:P'''''^Cp^^^^ClrK\\bgcd) + ^^Ol9{abgcd) ~ -^^CP''^'''Cp„''"Crnuvg{ab9cd) (B4c) 
5 

Vabcde = -^V(^abcd;e) (B4f) 

Here, we have used a slight modification of the basis of Riemann tensor polynomials suggested by Fulling, et al. 
[3^ in order to express the result in its simplest possible form. We have chosen an equivalent basis for Weyl tensor 
polynomials by ignoring all terms involving Rab and R and replacing Rabcd by Cabcd- We have also ignored all 
terms which vanish or are not independent under symmetrization of the free indices. Next, we have eliminated terms 
involving two covariant derivatives of a tensor in favor of terms involving two covariant derivatives of a scalar and 
single covariant derivatives of a tensor. We have done this as it is computationally faster and easier to calculate the 
covariant derivative of a scalar than of a tensor. Finally, we have used the identity [soj : 

V(r)C""^^(^V(s)C|pcrr|6) = ^.9(afc) C'"^^" V(5) CparK (B5) 

(where V(r) indicates r covariant derivatives) to combine some of the remaining terms. 

The motivation for working with V{x,x') rather than the Vn{x^x') is immediately apparent from the vast simplifi- 
cation that occurs. 



APPENDIX C: CALCULATING COORDINATE HADAMARD/DEWITT COEFFICIENTS TO 7"" 

ORDER IN THE SEPARATION OF THE POINTS 

In order to obtain expressions for the quasi-loca l eq uations of motion up to 7th order in the separation of the points, 
it is first necessary to extend the results of Ref. [3l[ by one order. As is noted by Ref. [lB| and previously by Rcf. 
[3^ . this is a relatively easy calculation due to the symmetry in the Green's function: 

v{x,x') =v{x',x) (CI) 

We use expression (|6.2p for v (x, a;') in the above to give: 

00 00 

v,,u (r) {t - t'f (cos (7) - 1)^' (r - r'f ^ (^') " (V) " 1)' " (C2) 

i^.k—Q i,j,k—0 

Now, in a similar vein to Appendix |^ if we take a sufficient number of symmetrized derivatives of both sides and 
then take the coincidence limit, it is possible to express odd order coefficients in terms of derivatives of the lower order 
coefficients. 
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For the cases of relevance to the current work, the calculation is made even simpler by spherical symmetry. The 
coefRcients Vijk are functions of r and M only (i.e. they have no dependence on t or 7). Additionally, since (t — i')^* 
and (cos(7) — !)■' are both even functions, they arc invariant under interchange of primed and unprimed coordinates. 
This means that we can express (jC2p in the simpler form: 

00 00 

E (0 - ^'f = E ^■'^■fc (^') - (C3) 

fc=0 k=0 

Taking r-derivatives of both sides of this equation and taking the coincidence limit r' ^ r gives the result 

Vijl = -^WyO.r (C4) 

For the present work, only the coefficients W301, V211, V121 and wqsi are needed. They are given here for convenience: 

1 AP {r - 2Mf (598rA/2 _ igsMr^ + 2Qr^ - SSSM^) 
-301 = (C5a) 

1 A/2 (r - 2M) (3352rA/2 - lOmMP + ll2r^ - 3240M3) 
- ^ ^15 (C5b) 

9 A/2 (228rA/2 - giA/r^ + 12r3 - 189A/3) 
„^2^ = V i C5 

1120 ^ 

1 A/2(-155A/r + 56r2 + 36A/2) 

3360 r' 
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